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Abstract 

We present an alternative way to calculate the screening of the static potential between two 
charges in (non)abelian gauge theories at high temperatures. Instead of a loop expansion of a gauge 
boson self-energy, we evaluate the energy shift of the vacuum to order e 2 after applying an external 
static magnetic held and extract a temperature- and momentum-dependent dielectric permittivity. 
The Hard Thermal Loop (HTL) gluon and photon Debye masses are recovered from the lowest lying 
Landau levels of the perturbed vacuum. In QED, the complete calculation exhibits an interesting 
cancellation of terms, resulting in a logarithmic running a(T). In QCD, a Landau pole in a s arises in 
the infrared from the sign of the gluon contribution, as in more sophisticated thermal renormalization 
group calculations. 



1 Introduction 

In quantum field theory, fluctuations of the vacuum give rise to the production of pair quanta which 
tend to screen (or antiscreen) the charge of a heavy test particle. If one perturbatively calculates the 
non-relativistic potential V(r) between two unlike static charges, say, in QED, the usual Coulomb-like 
behaviour is modified by the photon self-energy II (K 2 ) such that 

V(r) = [ d k e i% ' F — fl) 

[ ' J (2tt)3 k 2 + U(K 2 = -k 2 ) K ' 

where k — \k\ and K = (k°,k). Inserting the text-book result for H(K 2 ) and expanding for small 
distances k 2 3> m 2 , the quantum fluctuations lead to an effective coupling constant 



where A = exp(5/3)m e is a scale related to the electron mass m e . This is, of course, the familiar result 
of the running coupling in QED which is commonly obtained using renormalization group methods. In 
refs.[jl], pi, it has been shown that the running of a coupling constant at T = can be understood in 
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physical terms by the polarizability of the vacuum. The effects of fluctuations can be incorporated to a 
certain extent in a scale dependent dielectric permittivity e that defines an effective charge a e g = a/e. 
In vacuum, Lorentz invariance dictates that 

(ie = 1, (3) 

where /i is the magnetic permeability. Calculating fj.(k) at the momentum scale k and extracting the 
leading log contribution, one finally recovers the familiar expressions for the running couplings in QED 
and QCD. Then, asymptotic freedom can be interpreted in terms of a paramagnetic ground state. 



In this work, we extend the approach of to finite temperature and calculate an effective coupling 
constant a e s{k,T). Instead of a loop expansion, we evaluate the energy shift of the vacuum to order 
e 2 after applying an external (chromo)magnetic field H. The connection of magnetic permeability and 
dielectric permittivity at finite temperature is made by invoking a renormalization group argument. 
QCD with a magnetic background field at finite temperature has been studied in a number of works JjJ. 
In contrast to previous approaches, we lay out a non-technical calculation of charge screening without 
reference to propagators or self-energies, resorting to entities that have an immediate physical interpre- 
tation (energy densities and susceptibilities). Our work then allows an alternative, though slightly more 
phcnomcnological, view on screening at high temperature. 



A possible dissolution of bound quarkonia states, e.g. J /i/ 1 , was proposed long ago as an experimental 
signature of the quark-gluon plasma in heavy-ion collisions. The T-dependence of the interquark po- 
tential in QCD is therefore of great interest, and simulations of the potential in lattice gauge theory do 
indeed show a strong screening 0]. In perturbation theory, the quantity that enters the Fourier trans- 
form of the potential at finite temperature is the static limit of the longitudinal gauge boson self-energy 

Equivalently, one can define a dielectric permittivity by j(| 

„»,„_, + 3sSW2. © 

The perturbative one- loop thermal contribution to has been calculated long ago as JtJ: 

U L (0,k.T) = lJ- = {m e D f for QED, (6) 



and 



Ul ( 0. A • T) =[N c + ^A 9 -^-= (m c D f for QCD, (7) 



which defines screening masses mjj. Here, e and g are the electromagnetic and strong coupling, respec- 
tively, N c is the number of colours and Nf the number of thermally active flavours. Since the static 
limits of the self-energies are momentum-independent, the poles of the expression in (Q) are simply the 
gauge invariant Debye masses tud defined in eqs.(^) and (Q) and lead to an exponential damping of 
the potential V(r) ~ exp(— mr>r)/r. In particular, this form of 11^ has the consequence that gluons 
screen the strong interaction, in contrast to the zero temperature case, over long distances. However, 
the formula for the running QCD coupling constant commonly used in finite temperature calculations 
assumes that typical momentum transfers are of the order of the temperature, hence 

tts(T) ~ (UN c -2N f ) log(T)' (8) 
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In this expression, gluons therefore retain their antiscreening property, reflecting the ultraviolet sector 
of the theory. The transition to Debye screening is not obvious. Another troublesome feature of QCD 
screening at finite temperature is the behaviour of the Debye mass at next-to-leading order, which reads 

I 

mg ) (T) = m« + ^l Q g(£) ) (9) 

with rn^ given by eq.(Q). Here C is a constant arising from the ad hoc removal of infrared singularities 
involving chromomagnetic static modes. The appearance of the non-perturbative logarithmic term 
questions the applicability of loop calculations somewhat. Furthermore, whereas in QED the self-energy 
tensor is gauge independent, this is not the case in QCD, which makes the very definition of a Debye 
mass conceptually difficult. Finally, due to the nonlinear coupling of the gluons, relation (@) remains 
valid only within certain gauges (like temporal axial gauge) . An evaluation of the effective charge 
and its possible screening not relying on a Feynman graph expansion is therefore desirable. 

2 The zero temperature case 

In this section, we define our notation and briefly review the calculation of refs.Jl], g|. To obtain a 
scale-dependent fi, let us look at the change in the energy E of the vacuum when an external magnetic 
field H is applied: 

AE = ~[4ir X (H)]VH 2 ~E vac , (10) 

where fi(H) = 1 + 4ttx(H), and x{H) is the field-dependent magnetic susceptibility. As soon as the 
energies at H = and finite H are known to some approximation, a field-, or equivalently, scale- 
dependent fJ,(H) can be extracted. Later, the external field 2eH is identified with the scale K 2 at which 
the physical process is probed. 

For charged scalar fields, the general expression for the energy spectrum of a single Fourier mode reads 

= + 0. (n) 

distinguishing between particles (+) and antiparticles (— ). The dispersion relation LOk — k follows from 
the positive energy solution of the Klein-Gordon equation for massless, non-interacting particles. At 
T = 0, the occupation number nj: for the ground state is zero. Summing over particle and antiparticle 
states, we recover the familiar divergent zero-point vacuum energy Eq = Ylk Wk - F° r massless spin-i 
fermions, the energy without an external field becomes 

Ef = -2j2"k- (12) 

k 

The factor 2 arises from the spin summation, the factor —1 stems from the anti-commutation relation 
fermionic annihilation and creation operators obey. In the presence of the magnetic field H, we substitute 
— ► = dfj_ — igqAft, where q is the charge of the (anti)particle in units of the coupling g. Choosing 
the orientation of the i/-field along the z-axis, we construct a vector potential as — (Q,0,xiH,0). 
This choice for obeys d^A^ = 0. In the following, we treat QED and QCD in parallel and define 
e = qg. We have to solve for the energy spectrum of Uj)ijj{x) — 0, which is basically a relativistic version 
of the Landau theory for the diamagnetic properties of an electron gas. The solution for the energy of 
a single Fourier mode becomes 

E n ,k s , S3 = \Jki + 2eH (n + l/2 + s 3 ). (13) 
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In addition, the x\ space variable is shifted by —kz/ieH). Note that the energy depends only on two 
quantum numbers. The third is 'hidden' in the mentioned x\ shift. Here S3 = ±|, the ^-component of 
the spin. The Hs 3 term clearly shows the coupling of the spin to the external field, and hence, if the spin 
of the fermion is anti-parallel to the iT-field, the energy is lowered. For QCD, there is also an implicit sum 
over the colour charges q hidden in e = gq. Finally, for a vector gauge boson the H -independent energy 
is the same as for a scalar field, except that there is an additional factor of 2 counting the transverse 
spin degrees of freedom: 

E 3 =2^u k . (14) 

k 

The sum over colour degrees of freedom yields an additional multiplicative factor of N% — 1. In presence 
of the magnetic field, we separate the field into the classical background part A^ and the fluctuating 
quantum part A^. The equations of motion become D^G^ = 0, where G^ v is the gluon field strength 
tensor. With a suitable choice of background gauge, the energy for the two physical degrees of freedom 
of A^ can be written as 



En,k 3 ,s 3 = y/kl + 2eH(n +1/2 + s 3 ), (15) 

the same as in the fermionic case, but now with s 3 = ±1. Again, summation over the colour charges is 
implicitly assumed. 

We want to extract the leading log(-ff) contribution to the energy shift induced by the external field. 
With the total spin s of the particle considered and i — f,g: 

fel|fe2,fc3)«3 




K,k 3 ,s 3 = ^H + 2eH(n+l/2 + s 3 ). (17) 

Introducing a quantization volume V — L 3 , we replace the sum over ki and ^3 by an integral weighted 
with the density of states. Taking into account that the x\ variable was shifted, hi is restricted to 
< fc 2 < LeH. Then, 

y^^J dh Y, {eH - L) = ^ {eH) I dh - (18) 

To regularize the divergence, we will use a UV cut-off A such that < n < -^jj = «a and fc| < A 2 . The 
first idea would be to replace the sum over n by an integral. However, if we perform the shift n' = 2eHn, 
we find that the integral would be independent of H to leading order. That is, we would have recovered 
the vacuum result, in the absence of the field H. So what we need is the correction to the replacement 
of a sum with an integral. Such a correction term suitable for our case here is provided by a specific 
Euler's sum rule 



"2 nTl2 1 

]T /(n + l/2)= / f( x )dx--f'(x)r 



ni ■ 



We may now re-define the energy shift as 

~i + s 3 ) 



where 



{"A „,l A 
V f(n + 1/2 + s 3 ) - / dn f(n + 1/2 - 

fix) - ^(eH)£dk 3y /k% + 2eHx. 
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Since we are not interested in the soft modes of the order of eH (the leading-log behaviour is dominated 
by the UV behaviour of the theory), we split the sum into two pieces (TV <C ua) 

"A N n A 

E = E + E- 

n=0 n=0 n=N 

Let us treat S3 formally as a continuous variable. Taylor expanding in S3 (since n > N 3> S3), we are 
left with 

"A 



AE i = ( _ 1)2s J2 £ ( f( n + 1/2) - a 3 f'(n + 1/2) + f f"(n + 1/2) + 

s 3 n=N 

+ <f>{eH,N). 



(19) 



Now &(eH,N), which represents the contributions from soft modes only, does not depend on A. It is 
thus proportional to {eH) 2 for dimensional reasons, a small non-leading log contribution, and may be 
safely neglected. The linear term in S3 vanishes upon summation, and re-substituting e = gq, we find 



AE i = - 1 -V{gHf 



2tt 2 



1 



A 2 \ 



(20) 



The sum over a SU(N C ) multiplet of the squared charges q 2 is Nf/2 for the fundamental representation 
(Nf quark flavours) and N c /2 for the adjoint representation (the gluons). For QCD, the susceptibility 
becomes 

- // 

4ttx -> -g 



9 UN c -2N f , 

log 



A 2 



48tt 2 

which reproduces the expression obtained by renormalization group calculations if we identify 2eH = K 2 
For QED, the sum over the charge(s) is simply 1, so we obtain 

» 2 (2eH 



4ttx 



12tt : 



log 



V A 2 



again in accordance with eq.(g). Having outlined the calculation of 
part of the paper and switch on temperature. 



we now proceed to the main 



3 The temperature-dependent part 

At finite temperature T, the occupation number appearing in cq.(pT|) does not vanish anymore for 
the thermal ground state, instead nk — (exp(/3o;fc) — 1) _1 = tibe, the usual Bose-Einstein distribution 
function (/3 = 1/T). For fermions, = (exp(/3u>k) + 1) _1 = nFDi the Fermi-Dirac distribution function. 
Thus, when summing over the infinitely many degrees of freedom, we find for the total vacuum energy 
of a charged scalar field 



E° = }_^LU k (l + 2n BE (uj k )) 

k 



(21) 



The result clearly separates into the divergent vacuum part already treated and a finite, T-dependent 
part. In the case of a finite magnetic field H, the higher energy modes (17) are occupied with their 
respective thermal probabilities, and we can write (i — /,<?): 



AE 1 = 



E' 



K = 



E l -El 



E 

n,k 2 ,k 3 ,s 3 



E 



n,k 3 ,s 3 



(-iy s + 



k\ ,&2 5^3 -S3 



(-l) 2s + 



ex P (^, fc2 , S3 )-(- 1 ) 
2 



2.s 



and 



exp(/3w fc ) - (-1 



|2.s 



(22) 
(23) 

(24) 
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Again, we need to extract the leading thermal contribution to n(k). However, at finite temperature, 
relation (|^) does not hold any more. One could imagine to start with 

fi(k)e(k) =n(k) 2 , (25) 

where n is the index of refraction. This quantity is related to the photon or gluon phase velocity by 
v p — 1/n, and v p could be extracted from the (full) dispersion relation of the corresponding gauge boson 
since v p = ujk/k- However, eq.(|25|) holds only for "on-shell" propagating gauge bosons. Since Lorentz 
invariance is formally broken by the presence of the heat bath, /i and e become functions of fc° and fc, 
and eq.(|2^) reads, more explicitly, fi(oJk, k)e(ujk, k) = n{k) 2 . Here, we consider an off-shell external field, 
so a relation between /i(0, k) and e(0, k) is required. 



The total energy density of the system can be written as the sum of the field and the induced medium 
energy density: 

£tot = ^ 2 + ^$>£\ (26) 



2 V 

Introducing an effective field H c s, we rewrite £ as 

, _ 1 2 _HeH) 
c-tot — izti. 



tot - g -off - 2~i*p- 

In the last step we made use of the fact that eH has to be a renormalization group invariant, so 
eH = e e ffiJ e ff. The effective coupling constant is now defined by §, || 

using ( |l0| ) and (]2^). Replacing 2eH by k 2 , as at T — 0, our master formula hence reads 

a cS {k,T) = — — = a (27) 

e(k,T) l-47rx(fc,T) 

The thermal piece of eq. (|23|) can be compactly re- written as 



(oo °° 
Aee f dx 



sjx 2 + b(n+ 1/2 + s 3 ) 



S3 n=0 ' 



exp ^x 2 + 6(n + 1/2 + s 3 ) - (-1 



(28) 



where x is dimensionless and b = 2eH/T 2 is a measure for the ratio of quantum and thermal effects. We 
consider the high-temperature limit b <C 1 in the rest of the paper. 



3.1 A first (incomplete) approximation 

The sum appearing in expression ( p8| ) obviously cannot be evaluated exactly. It is instructive to work 
out the first intuitive approximation to the sum although we will show in the next section that it is too 
crude. 



Consider the fermionic part. Note that the factor b(n + 1/2 + S3) plays the role of a mass term in the 
integral in eq. (p8|) , so the contribution of the terms in the sum becomes exponentially suppressed as 
n increases. In contrast to the T = case we are therefore interested in the behaviour of the sum at 
small n where the S3 spin component is not negligible. Thus we cannot apply a Taylor expansion in S3, 



G 



as done in (JTq|) , and need an exact summation over S3. Isolating then the lowest lying Landau mode 
(n = 0, S3 = —1/2) and combining the remaining expressions into a single sum, we find 



VT 4 



E f = l + ±y 

24 i 2 ^ 

n=0 



dx- 



y/x 2 +b(n+l) 



exp ( y/x 2 + b{n + l) \ + 1 



(29) 



Since b <C 1, the terms in the sum vary slowly with n, so we can again try to trade the sum for an 
integral over n: 



~ f b 2 

E = — — 4 



24 



dr 



r 2 \/r 2 4- b 



7T 2 J exp (s/r 2 + b) + 1 ' 



Neglecting terms of order b in the integral, we obtain as a first approximation 



E: 



1, 



(0) 



-VH 2 



(rn 



e \2 

d) 



2eH 



7 7T 2 

4 15 



VT 



(30) 



(31) 



for QED with the HTL Debye mass defined in eq.(^). The second term is simply the energy Eq(T) of 
a thermally excited, non-interacting massless fermion-antifcrmion pair, i.e. the thermal energy of the 
unperturbed vacuum that has to be subtracted, cf. eq. (p2|) . This means that we have recovered within 
our simple framework the perturbative one-loop HTL result from the lowest Landau level contribution 
to the energy of the magnetically perturbed thermal vacuum. The energy difference that enters in ( |To| ) 
already yields Attx{H, T) as the expression in square brackets, and the effective coupling constant reads, 
following eq.(|27|), 

a eS (k,T) = 



a 



e(k,T) 



(mh) 2 
k 2 



as within HTL perturbation theory, eq.(|^). 



For QCD with Nj flavours, we obtain 



E 



(0) 



--VB* 



" L DJ 

2eH 



with the fermionic part of the squared QCD Debye mass (Q) , 



7tt 2 , 
Nf W VT 



(32) 



Nf/6 g 2 T 2 . For the total evaluation 



of the QCD susceptibility, we need to add the contribution from the gauge bosons. At zero temperature, 
contributions from "unphysical" gluon states in the calculation of the energy spectrum, eq.(|l5|), are 
exactly cancelled by Fadeev-Popov ghost contributions within the background gauge condition used 
here. Since we only consider excitations of energy levels that were evaluated at T = 0, no ambiguity 
arises and we still work only with physical gluon degrees of freedom with two polarization states. We 
proceed in close analogy to the fermionic case: first, we sum over S3 = ±1. A subtlety arises since the 
combination n = and S3 = — 1 in cq.(p8|) issues a negative value under the square root for small x. 
This 'tachyonic' mode is related to a possible instability of the vacuum [^0). Its effect on the magnetic 
field over large distances, however, is negligible, the sum over n for S3 = — 1 therefore starts only at 
n = 1. Isolating again the lowest lying physical Landau level (n — 1, S3 = —1) contribution to the sum, 
we are left with 



E 9 
VT 1 



E 9 



b 
2^2 



dx- 



exp 



6/2 



„=o I exp ( y/x 2 + b(n + 3/2) 



(33) 
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Replacing the sum by an integration, setting b = in the integrals and summing over colour, the result 
becomes 



R 9 
^(0) 



1, 



-VH 2 



2eH 



+ 2(N 2 -1)^-VT 4 
15 



(34) 



Again, the last term is the thermal energy Eq(T) of the unperturbed SU(N C ) gluon vacuum. The 
expression in square brackets exactly corresponds to the gluonic part of the squared QCD Debye mass, 



l D,g 



N c /3 g T . Putting all pieces together, the effective coupling becomes 



a s ,es(k, T) 



1 



(m c D ) 2 ' 
k 2 



very similar to the QED case. In our model, the Hard Thermal Loop (chromo)electric Debye masses 
therefore appear as the lowest Landau level contribution to the energy difference that arises when one 
probes the thermal vacuum by a (chromo) magnetic field. It is worth noting that, in this approximation, 
the alignment of an external field always increases the thermal energy of the vacuum, regardless of the 
non-abelian structure of the theory. Therefore xi^i E) is always negative and we conclude, using eq.(27), 
that the static potential would become screened by both fermions and by gauge bosons. 



3.2 A better approximation 

However, additional contributions to eqs. @, © and [M of the same order in e 2 arise from two 
sources. First, the expansion of the integrals (30|) and (33) in b is similar to the high temperature 
expansion of loop integrals with massive particles in the small quantity mo/T. The appendix contains 
the relevant formulae. Second, the correction to the replacement of the sum by an integral yields terms 
to order b and b 2 that are provided by the Euler-MacLaurin formula 



N " 
n=0 I 



f(x)dx + - [f(N) + /(0)] + - [/'(AO - /'(0)] 



(35) 



where the dots denote terms with higher derivatives in f(n). For our purposes, eq. (|35|) , taking A —> oo, is 
sufficient, as long as f(x) £ C 2 for x £ [0, A]. When calculating the thermal contribution to the vacuum 
energy, we include these correction terms to the integral in the following and expand all integrals in the 
small parameter b, using the relations presented in the appendix. The summation of all terms to order 



e then alters the results in (B2h and 



qualitatively. 



3.3 Results for QED 

For fermions, we start with eq.py). Defining 5 2 — b, we obtain to order 5 4 



E 



(i) 



r 2 \Jr 2 + S 2 



S 2 2 . , 

1 / dr 

24 7T 2 ./ exp(vV 2 + 5 2 ) + 1 ' 2 



1 

12 



o 
2ir 2 



dx 



Vx 2 + 5 2 1 + cxp(-Vx 2 + S 2 ) J exp(Vx 2 + S 2 ) + 1 



S 2 



~* J exp(Vr 2 + 5 2 ) + 1 



(36) 



Using the functions /, (y) and g + (y) defined in the appendix, we re- write 

S 2 2 . _ 5S 2 . /rt S 4 „ ,„ S 4 



E 



f _ _ 

(i) 24 



2tt 



24tt 2 



(37) 
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Expanding in S and keeping all terms up to 0(S 4: ), surprisingly all terms of order S 2 cancel, and we are 
left with 

%-^ + 5S» k *(^ < 38 » 

with 7 = 0.5772... and the constant Af — exp(27 — 1)/tt 2 ~ 0.12. The first term is the well-known 
thermal part of the vacuum energy in the absence of the field H . Since S 2 -C 1, the alignment of a 
magnetic field hence decreases the energy of the vacuum at finite temperature, in contrast to the result 
of the previous section. The susceptibility in QED therefore becomes 

A-K X {H,T) = -^-\og{A f ^j. (39) 

Note that the pre-factor of the logarithm is the same as in the zero-temperature case! Including this 
pure quantum correction for the QED running coupling, eq.(|J), all momentum dependence drops out, 
and we finally obtain with ( p7| ) 

a e 7+7/6 a 

a cS (k,T) = with A T = m e ~ — , (40) 

1 a ( T z \ 7T 3 



which is valid to order e 2 and for momenta 



m e < k < T. (41) 



Using the zero-temperature coupling a c s(k) from eq.(|2]), the effective coupling can, to this order, be 
rewritten as 

a cS (k, T) = a ef f (T) - a cS ((k) ~ 3 T). (42) 

So the common practice used in perturbation theory to simply take the running of the coupling at zero 
temperature and set as the scale the thermally averaged momentum scale (k) ~ 3 T does indeed find 
support from our calculation for QED. 



3.4 Results for QCD 

For QCD, the fermionic contribution takes a form similar to the QED result, 



T 2 



(43) 



Note that the pre-factor of the logarithm is again the same as at zero temperature. The calculation of 
the gluonic part of x runs along the same lines. Starting with eq.(^) and setting S 2 = 6/2, we obtain 
by use of the functions hi(y) and g~{y) 



E 9 = ^h 5 (V3S) + ^ \7h 3 (V36) + h 3 (5)] + ^ 



hi{S) + 3hi(y/3g)-^g-(Va5) 



Expanding, we find that all terms of order 5 3 cancel and the result becomes 



E 9 = 



~15 



b 

12 



5b 2 



128tt ; 



■ log (A g b) , 



(44) 



(45) 



with the constant A g = exp(27 — 13/10 + 11/5 log3)/(327r 2 ) ~ 0.03. Similar to the fermionic part, the 
alignment of a chromomagnctic field hence always lowers the energy of the gluonic vacuum, however, 
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Figure 1: The ratio a StC s/a s of eq.([47|) as a function of k/T for a coupling g = 0.1 and Nf = (solid 
line). For comparison, the corresponding ratio in the HTL calculation, using the Debye mass of eq. (]?]), 
is also plotted (dashed line). 



this time the difference goes linearly in b and not only logarithmically. Finally summing over colour, the 
gluonic susceptibility reads 



jv c .9 2 r 2 

3 2eH 



5N C 
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'48tt 



log A 



2eH\ 



In this expression the sign is reversed as compared to eq.(f43|), and 4ttx is always positive, 
thermal result for QCD, excluding the T = contribution, becomes 



(46) 



The total 



ot s ,es(k, T) 



a.s 



47riV c T z 
'~3~~k* 



12tt 



5AUoe 



A g k 2 

rp2 



2N f log 



T 2 



(47) 



Since the gluon contribution dominates by far over the logarithmic fermionic, there is antiscreening at 
high temperature and long distances. This result is in contrast to expectation and lattice results on the 
interquark potential Extrapolating eq.(^7|) beyond the kinematical region g <C k/T 1 where our 
approximations are valid, a Landau pole appears in the infrared region k/T ~ g. Figure 1 shows the ratio 
Q<s,eff/ 1 ot s as a function of k/T for a weak coupling g = 0.1, compared to the usual HTL result. A similar 
behaviour is also found in more sophisticated renormalization group analyses of the running coupling 

p"2| , |l3]|). We note that our results compare quite well with the 
|12| . Since eq.([47|) depends only on the dimensionless quantity k/T, 
taking the limit T — > oo at large k is in a sense equivalent to probing the infrared region k — > at smaller 
T, indicating that non-perturbative physics plays an important role even at high T. The necessity for a 
soft magnetic mass ~ g 2 T as an infrared regulator in loop calculations (see eq.(^|)) supports this line of 
reasoning. 



at finite temperature (see, e.g. 
numerical solutions obtained in ref. 
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4 Conclusions 



We have presented an alternative way to calculate the screening of the static potential between two 
charges in (non)abelian gauge theories at finite temperature by looking at the magnetic properties of 
the vacuum. Instead of a loop expansion, we have calculated the energy shift of the vacuum at finite 
temperature to order e 2 after applying an external (chromo)magnetic field if as a probe. Magnetic 
permeability and dielectric permittivity have been connected by a renormalization group argument. 

Using a high-temperature expansion H/T 2 <C 1, the gluon and photon Debye masses appearing in 
the HTL calculation have been recovered in a first, though incomplete approximation, originating from 
the lowest lying Landau level contribution to the thermal energy. Taking into account all contributions 
to order e 2 in QED, the final expression in the kinematic region ( [4l|) shows a logarithmic, momentum- 
independent running of a with temperature, as expected from simply inserting the average thermal 
momentum (k) ~ 3 T in the zero-temperature running coupling. 

In QCD, we have found indications for a Landau pole at small k/T that arises, as in more sophisti- 
cated thermal renormalization group calculations, from the sign of the gluon contribution, despite well 
controlled approximations and a completely different approach as compared to conventional perturbation 
theory. Our calculation may serve as yet another indication that an expansion in a presumably small 
coupling g at high temperatures ceases to yield sensible results for some quantities, and that this failure 
is not specific to a Feynman graph expansion. Truly non-perturbative input that is probably linked to 
the understanding of confinement is then called for. 



5 Appendix 



Here we present the formulas used to evaluate the (S-dependent integrals appearing in eqs.(44) and (|3j 
Fermionic integrals of the form 



DC 

'n(y) = J dx 





\J x 2 + y 2 exp(^/ x 2 + y 2 ) + 1 
can be expanded for small y (note our slightly different convention compared to |))). Especially, 



My) = 4[ lo §(f) ' 



My) 



tt 2 y 2 
12 + ~A 
7tt 4 tt 2 



logg)+ 7 -i 



and 



7 "4 



where 7 — 0.5772... is the Euler-Mascheroni constant. For bosons, 

„n— 1 1 



K{y) 



The corresponding expansions read 



dx 



x 







\Jx 2 + y 2 exp(y/ x 2 + y 2 ) — 1 



i/1 7T 1 

- Yy + 2 



, , , TT' TT y 

r 2 



V47T 
2 



log 



\AttJ 



hh(y) = v H — v H y 



•7"2 



Hi) 



and 



■7 

' 4 



+ . . 



(48) 
(49) 

(50) 



(51) 
(52) 
(53) 
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For the evaluation of derivative terms, we need the leading log(y) behaviour of integrals such as 



g ± (y) = gf(y) + gf(y) 



dx 



x 2 + y 2 1 ± exp(— x 2 + y 2 ) J exp(^/a; 2 + y 2 ) ± 1 



(54) 



The expansion of the first term in brackets, g 1 , is known since g^ — f\ and g 1 = h%. For the evaluation 
of the second a trick is convenient. Introduce a parameter a to write 



gf (y; 



dx- 



exp(a-\/ x 2 + y 2 ) 
(exp(a^x 2 + y 2 ) ± l) 2 



(55) 



Obviously, gf(y; 1) is the sought quantity. Now gf(y; a) can also be written as 



d_ 

da 



dx 



\/ x 2 + y 2 exp(a-\/ x 2 + y 2 ) ± 1 
Expanding eq.(^35[) for small y hence yields 

g^iv, = 20 — ^ o ^ ' • ' f° r bosons and 

g^{y] cx) = —- h • • ■ for fermions 



1 

2a 



Setting a = 1 and putting the pieces together, the leading-log behaviour of eq.(p4|) is 



g + (y) = -o[log(- )+7+i 



+ 7 + 1 



2 



for fermions and 
for bosons. 



(56) 

(57) 
(58) 

(59) 
(60) 
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